Bloch-like Electronic Wave Functions in Two-Dimensional Quasicrystals 
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Electrons in quasicrystals generically possess critical wave functions that are neither exponentially- 
localized nor extended, but rather decay algebraically in space. Nevertheless, motivated by recent 
calculations on the square and cubic Fibonacci quasicrystals, we investigate whether it is possible to 
obtain extended wave functions expressed as linear combinations of degenerate, or nearly-degenerate, 
critical eigenfunctions. We find that not only is this possible, but that the wave functions that 
emerge are Bloch-like, exhibiting the quasiperiodic long-range order of the underlying quasicrystal. 
We discuss the significance of this result for the study of electronic properties of real quasicrystals. 
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Although research in quasicrystals celebrated its 
twenty-fifth anniversary last year [l|, Q , many gaps still 
remain in the fundamental understanding of their elec- 
tronic properties A major shortcoming is the lack of 
a quasiperiodic analog of the Bloch theorem, which plays 
an important role in the study of the electronic proper- 
ties of periodic crystals, serving as the basis for band- 
structure and electronic-transport theories. Naively, one 
would expect to encounter a quasiperiodic version of the 
Bloch theorem when dealing with quasicrystals [1, 
predicting the existence of extended eigenfunctions that 
display the underlying quasiperiodic long-range order. 
Some experiments even seem to indicate that such ex- 
tended eigenfunctions exist in real quasicrystals [6] . How- 
ever, as previous studies have shown 0, Isl9]> the 
electronic eigenfunctions in quasicrystals are generically 
found to be critical, with algebraic spatial decay. 

We have recently reintroduced the use of separable 
models — the square and cubic Fibonacci tilings [Io|— 
for studying the physical properties of quasicrystals (for 
precise definitions of the terms 'crystal' and 'quasicrys- 
tal' see [lH and [II]). These models allow us to use 
well-known results for the 1-dimensional (Id) Fibonacci 
quasicrystal and extend them to higher dimensions. Our 
past studies of the electronic energy spectra of these qua- 
sicrystals [l3|, EH, [HI led us to conjecture (as we explain 
below) that linear combinations of degenerate eigenfunc- 
tions, all of which are critical, may produce extended 
wave-functions, but we failed to find such combinations. 
In this Letter, using a new iterative search algorithm for- 
mulated by Elser et al. [T^j, we establish that not only 
do such wave-functions exist, they actually exhibit Bloch- 
like properties. 

The square Fibonacci tiling is constructed by super- 
imposing two Id Fibonacci grids perpendicular to each 
other, as shown in Fig. [T] The off-diagonal tight-binding 
model on this quasicrystal 
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on-site energies (taken to be zero), and hopping that 
is restricted along tile edges, with amplitude 1 for long 
edges and T for short edges, where we take T > 1. The 



FIG. 1: A section of the square Fibonacci tiling, whose geo- 
metric properties are described in Ref. [l(J. 



Schrodinger equation is then given by 

T i+1 ^>(i + lJ)+T^(i-l,j) 

+T j+1 H>(i,j + 1) + Tj^(i,j - 1) = £*(*, j), (1) 

where ^(i,j) is the value of a 2c? eigenfunction on a ver- 
tex labeled by the two integers i and j, and e is the 
corresponding eigenvalue. The hopping amplitudes Tj 
are equal to 1 or T according to the Fibonacci sequence 
{Tj} = {l,r,l,l,T,l,T,l,l,T,l,l,r,l,r,...}. The 
fact that Eq. |T]) is separable allows one to use the known 
solutions for the Id problem 0, @, ©, \t& 0, 0, EH in 
order to study higher-dimensional models (as was done 
for similar models in the past [12, El El EE ES El El ) • 
Two-dimensional eigenfunctions are therefore Cartesian 
products of Id eigenfunctions, *& n ,m = ipn ® VVn, and the 
corresponding 2d eigenvalues are the pairwise sums of Id 
eigenvalues e n m = e n + e m . The Id eigenvalue problem 
is solved by considering periodic approximants of the Id 
Fibonacci quasicrystal. The unit cell of the N th order 
approximant consists of Fn sites — where F/v is the N th 
Fibonacci number, defined by Fjsr — Fj\r—i + -F/V— 2 with 
Fq = F\ = 1 — and therefore its spectrum consists of Fn 
bands. Explicit diagonalization of the Hamiltonian for a 
single Id unit cell, with periodic or antiperiodic bound- 
ary conditions, yields the edges of these bands, and hence 
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provides the full structure of the Id spectrum. 

Our studies of the energy spectrum of the 2d model 
revealed that for T close to 1 (i.e. near the periodic 
limit) the addition of Id spectra gives rise to a 2d spec- 
trum containing energy intervals, similar in its proper- 
ties to the spectrum of a periodic crystal, even though 
the Id spectrum is nowhere dense, containing no inter- 
val. This behavior of the spectrum led us to expect that 
eigenf unctions would be extended, even though they are 
constructed as Cartesian products of critical eigenfunc- 
tions. We proposed that the degeneracy of the 2d spec- 
trum, at values of T close to 1, may provide a mecha- 
nism for the emergence of extended wave functions as 
the dimensionality of the quasicrystal increases. Our 
conjecture was based on the idea that appropriate linear 
combinations of degenerate algebraically-decaying eigen- 
functions, centered about different sites but with suffi- 
cient overlap, may yield an extended function that spans 
the whole quasicrystal. To test this hypothesis we used 
an ad hoc approach to try to find an extended linear 
combination of zero-energy eigenfunctions [3] • The fact 
that the zero-energy eigenfunctions are macroscopically- 
degenerate, gave us great freedom, yet we only managed 
to obtain wave functions that were extended along one 
direction, but very strongly localized in the perpendicu- 
lar direction. 1 

The iterative search algorithm of Elser et al. [l|| gives 
us a new way of maximizing the extent of wave func- 
tions that are spanned by a given set of eigenfunctions. 
The algorithm is based on formulating a problem as an 
attempt to simultaneously satisfy two constraints, rep- 
resenting them geometrically, and using projections iter- 
atively between the representations until convergence is 
obtained. In our problem, the first constraint is for the 
function to be spanned by the given set {$> n .m\ of degen- 
erate eigenfunctions, which determines a vector subspace 
and the natural projection operator into it. The second 
constraint is for the wave function to be totally extended, 
with equal amplitudes on all lattice sites. This constraint 
determines a projection operator in which the amplitudes 
of a wave function are all set to be equal, while its phases 
are kept unchanged. The two operators used for the ap- 
plication of the algorithm are thus 
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(3) 



Fjv|*(j.*0I 

and the search is stopped after applying the projection 



1 Note that in Ref. [T3 | — where we showed that for every Id 
eigenfunction ip n (i) there exists a 2d eigenfunction ^> n {i^j) = 
(— 1)3 il) n (i)ip n (j) with zero energy — we gave an incorrect esti- 
mate of the spatial extent of fy n . 
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FIG. 2: The inverse participation ratio as a function of iter- 
ation steps for different values of the relative hopping ampli- 
tude T and different types of initial guesses. 



operator Pi , ensuring that the final wave function can be 
spanned by the required set of eigenfunctions, although 
generally it will not be totally-extended. 

The algorithm minimizes the £2 distance from the set of 
totally-extended wave functions. However, a more com- 
mon measure for the extent of a wave function is its in- 
verse participation ratio, 



/{*} = 



E 3 -, fc l*(j.fc)l 4 
(£ M l*(j.*0l a 



2 - 



(4) 



which we use here to trace the convergence of the algo- 
rithm. We find that the iterative search algorithm quite 
easily succeeds in finding extended linear combinations 
for all three kinds of degenerate sets considered below. 
Fig. [2] shows the inverse participation ratio as a function 
of iteration number for some typical examples. The al- 
gorithm is seen rapidly to converge to a steady, minimal, 
value of the inverse participation ratio, corresponding to 
a maximally-extended wave function. Naively, one may 
expect the convergence to accelerate if one picks the most 
extended eigenfunction as an initial guess, but it turns 
out that for low values of T the algorithm cannot escape 
the neighborhood of such initial guesses. In such cases, 
an extended function with random phases, taken as the 
initial guess, yields a better final outcome, as demon- 
strated in Fig. [21 

Turning to actual results, we start by reexamining the 
possibility of generating an extended wave function by 
using the macroscopically-degenerate zero-energy eigen- 
value. Applying the iterative search algorithm to the 
set of zero-energy eigenfunctions yields highly-extended 
eigenfunctions, despite our past failure [TJ|, as can be 
seen in the supplementary material [29| . 

The high degeneracy of zero-energy eigenfunctions is 
an artifact of our model which does not have any real 
physical significance. We therefore consider non-zero en- 
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FIG. 3: The largest number of bands with a common energy 
value divided by Fn as a function of the relative hopping 
parameter T, for different orders of approximants. The con- 
vergence of the curves implies that the spatial density of de- 
generate eigenfunctions, per unit area, decreases as the order 
of approximant increases. 



ergies e > 1 + T, taken from the upper quarter of the 
2d energy spectrum, which is fully-contained within the 
interval —2 — 2T to 2 + 2T. It is possible for certain 2d 
eigenvalues to be obtained as sums of pairs of Id eigenval- 
ues in a number of different ways, and thus have nontriv- 
ial degeneracy. To approximate this degeneracy for finite 
order approximants, we consider for each energy in the 
2c? spectrum the number of bands — constructed as sums 
of pairs of Id bands — in which it is contained. We then 
pick a representative eigenfunction from each of these 
bands to be included in the degenerate set. Applying 
the algorithm on such sets of degenerate eigenfunctions 
yields extended wave functions with some appearance of 
long-ran ge o rder, as can be seen in the supplementary 
material [29(. 

Despite this apparent success, one should be careful in 
assigning it any physical significance, because for a real 
electron to be able to explore the degenerate set of eigen- 
functions these functions must possess sufficient spatial 
overlap. Fig. [3] plots the number of bands that contain 
the most degenerate energy, e > 1 +T, as a function of T, 
divided by the corresponding Fibonacci number Fn for 
four consecutive orders of approximants. The steep de- 
cline of the curves as T increases shows that degeneracy 
is negligible for high values of T, as expected. The con- 
vergence of the curves for increasing N implies that the 
degeneracy is proportional to Fn whereas the number 
of sites is equal to F N . Thus, as the order of approxi- 
mant increases, the spatial density of degenerate eigen- 
functions decreases. In the quasiperiodic limit the degen- 
erate eigenfunctions might only have negligible overlap, 
and a real electron may retain the critical nature of its 
eigenfunction. 

Finally, we consider the most physically-relevant situa- 



tion in which thermal fluctuations, a small amount of dis- 
order, or the existence of inhomogeneous external fields, 
impose a small uncertainty in the exact energy of an elec- 
tron. We take this uncertainty to be Ae = 0.01, which 
is less than a quarter of one percent of the total band- 
width. We treat all the eigenfunctions with eigenvalues 
within a window of size Ae as if they were degenerate, 
and apply the iterative search algorithm to the largest 
nearly-degenerate set, with energies above 1 + T. The 
upper row of images in Fig. [4] show the intensities of the 
wave functions that are obtained by the algorithm, for 
three values of T. Even though the algorithm is devised 
only to maximize the extent of the wave functions, re- 
gardless of their actual structure, for sufficiently strong 
quasiperiodicity, or large T, it gives rise to the emergence 
of Bloch-like wave functions. 

For relatively low values of T the spatial decay of the 
critical eigenfunctions is relatively slow and the overlap 
between them is such that the wave functions obtained 
by the iterative algorithm appear to be smeared over the 
entire approximant with very weak apparent long-range 
order. As T increases the spatial decay of the critical 
eigenfunctions becomes more dominant, and the obtained 
wave functions display Bloch-like behavior, as can be seen 
in the plots for T — 2.5 and T = 5. These wave func- 
tions are Bloch-like in the sense that their intensities 
| <3> | 2 exhibit the underlying quasiperiodic order of the 
quasicrystal, although one cannot associate with them 
a well-defined non-reciprocal-lattice wave-vector, which 
is an additional requirement of the Bloch theorem (see 
Refs. [30] and [3l| for a discussion of reciprocal lattices 
of quasicrystals). 

This observation is supported by the Fourier trans- 
forms of the probability densities l^l 2 , shown at the bot- 
tom of Fig. [U Each Bragg peak with amplitude above 
some threshold is represented by a circle whose radius is 
proportional to the amplitude. While some Bragg peaks 
appear already for T = 1.4 they lack square symmetry, 
and the diffraction contains a considerable diffuse compo- 
nent (not shown due to the threshold) . For higher values 
of T it clearly shows Bragg peaks in a pattern similar to 
that of the square Fibonacci quasicrystal itself. 

To see whether these results are valid in the quasiperi- 
odic limit Fig. [5] plots the number of eigenfunctions in 
the most populated energy interval Ae as a function of 
T for approximants of order 9 — 12 (Fn = 55 — 233). 
The number of degenerate eigenfunctions is divided by 
the number of sites, F N , and can be seen to scale with 
it by the convergence of the different curves. This im- 
plies that the spatial density of nearly-degenerate eigen- 
functions remains constant, so that appropriate linear 
combinations can retain the Bloch-like nature as one ap- 
proaches the quasiperiodic limit. 

The existence of Bloch-like wave functions in the 
square Fibonacci quasicrystal encourages us to seek out 
similar results in more realistic models of quasicrystals, 



FIG. 4: The intensities |^| 2 of the wave-functions obtained by the search algorithm on a 144 x 144 approximant for different 
values of T are shown in the upper row. The underlying Fibonacci order becomes evident as T increases. This is supported 
by the emergence of Bragg peaks in the diffraction patterns, shown on the bottom row. Bragg peaks are represented by 
circles whose radius is proportional to their amplitude. Additional images and animations can be found in the supplementary 
material [29t ] . 



7 1 




1 I 1 1 1 1 1 1 1 L- 

1 1.5 2 2.5 3 3.5 4 4.5 5 



FIG. 5: The number of eigenfunctions in the most populated 
energy interval, in the upper half of the 2d spectrum, divided 
by the number of sites in the 2d quasicrystal (F^) as a func- 
tion of T. The convergence of the curves implies that the 
spatial density of nearly-degenerate eigenfunctions does not 
change as the order of approximant increases. 



like the Penrose tiling. Even though individual eigen- 
functions are critical in these models, we intend to inves- 
tigate whether combinations of nearly-degenerate eigen- 
functions can be extended. These nearly-degenerate com- 
binations of functions may turn out to play an impor- 
tant role in realistic situations, significantly affecting the 
transport properties and quantum dynamics of electronic 
wave packets in real quasicrystals. 
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